Analysis of Measured Magnetization 
and Pull Characteristics 

By R. L. PEEK, Jr. 

(M.iuusrript ret-eived September 21, 1953) 

It is shown in this article that the observed magnetization relations of 
most ordinary electromagnets conform to simple expressions which can be 
interpreted as the flux vs. magnetomotive force equations of a reluctance 
network, analogous to the current-voltage equations of a resistance network. 
To the extent of such conformity the magnetic circuit constants characterizing 
the network suffice for the evaluation of the field energy and pull charac- 
teristics of the electromagnet. The agreement between the observed mag- 
netization and these simple relations is close in the region of linear mag- 
netization, and is adequate for engineering purposes at higher flux densities, 
but the extent of agreement in the latter range varies with the type of struc- 
ture and the location of the magnetic parts which first approach saturation. 

Specific analytical and graphical procedures are given for the evaluation 
of the magnetic circuit constants from both pull and magnetization measure- 
ments. These procedures employ relations which give linear plots imlicating 
the degree of conformity of the observed relations to the expressions used to 
fit them. The relation of the measured constants to those which can be esti- 
mated in design is discussed, as is the use and application of the measured 
constants in development and engineering studies. 

1 INTRODUCTION 

III the design of telephone relays and similar smtching apparatus, the 
characteristics of the electromagnet which serves as motor element may 
be distinguished from those of the mechanical system of contact springs 
and actuaring meml>ers which it operates. The performance of the elec- 
tromagnet is characterized statically by the mechanical work done for a 
given coil energization, and dynamically by the time required to actuate 
the mechanical load, inclnding in this the hiertia of the moving parts. 

Both the potential work output of the electromagnet and the energy 
stored in developing its field can be evaluated from its magnetization 
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relations. The ronsequent dependence of the pull and timing characteris- 
tics upon the magnetization relations makes the measurement and analy- 
sis of the latter fundamental to the understanding of performance and 
its relation to design. 

Procedures are described in this article for the evaJuation from meas- 
ured magnetization relations of a few parameters which suffice for the 
determination of the pull and of the field energy under given conditions. 
These parameters, the magnetic circuit constants, characterize the elec- 
tromagnet to which they apply. They are used as measures of per- 
formance in comparing different structures, or in studying the effect of 
dimensional or material variations in a particular design. 

In addition to their significance as parameters summarizing measured 
magnetization relations, the magnetic circuit constants may be inter- 
preted as the observed values of quantities postulated in the magnetic 
circuit approximation to static field theory. This approximation is dis- 
cussed in a companion article,' which describes methods of estimating 
the values of the magnetic circuit constants from the configuration, 
dimensions, and materials constant of the electromagnet. Comparison of 
observed and estimated values of these constants has served as a guide 
in developing these methods of estimation. A complete design method- 
ology is provided by the ability to both estimate and measure the mag- 
netic circuit constants characterizing the performance of electromagnets. 

In the experimental evaluation of the magnetic circuit constants de- 
scribed in this article, the basic method is that employing measured 
magnetization relations. Because of the dependency of the pull upon the 
magnetization relations, pull measurements may also be used, subject to 
certain limitations, for the evaluation of the magnetic circuit constants. 
The article includes description of procedures for doing this. 

The notation used in this article conforms to the list that is given 
on page 257. 

2 MAGNETIZATION RELATIONS 

The magnetization relations of an electromagnet give the average 
flux <p linked per turn of the winding as a function of the two determining 
variables: appfied ampere turns NI and armature position x. The rela- 
tions are usually shown, as in Fig. 1, as a family of curves giving (p 
versus NI for various values of x. Armature motion is usually rotary, 
and the choice of the location at \\'hich x is measured is a matter of con- 
venience. The curves of Fig. 1 apply to the electromagnet shown in Fig. 
2, the AJ (heavy duty) tyi)e of wire spring relay described in an article 
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Fig. 1 — Magnetization curves. 

by A. C. Kellei".^ In this case it is convenient to measure x at the center 
line of the actuating card, whose location with respect to the axis of 
rotation is indicated in Fig. 2. 

Apart from interpretation in terms of magnetic theory, any observed 
point on a magnetization cur\'c represents a measurement of the elec- 
trical energy U stored in the coil, as given by the integral 



Jo 



ei dt, 



where c = d{N<p)/dt. It is represented by the area between the magneti- 
zation curve and the axis of ip, and is given by: 



U 



Jo 



NId<p. 



(1) 



This stored field energy may, in principle, be recovered in the decay 
of the field, except for the portion lost by hysteresis in the material, 
represented by the area between the increasing and decreasing mag- 
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Fig. 2 — Magnetic parts, coil, and actuating card of AJ relay. 

netization curves. The dashed curve of Fig. 1 is the decreasing magneti- 
zation curve for the same gap as the adjacent magnetization curve. 

[/ is a function of x, and changes as the armature moves. If this 
motion occurs at a constant value of ^, there is no induced voltage and 
hence no electrical transfer of energy. A decrease in U must therefore 
represent mechanical work done on the armature. For a differential 
change dx in armature position, the mechanical work is Fdx, where F is 
the pull exerted on the ai-mature. A general expression for F is therefore 
given by : 



F = 



dU 

ox (^i^onatant) 



(2) 
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It follows from these considerations of energy balance that the mag- 
netization relations determine the relations among the electrical input 
to the coil, the stored energy, and the mechanical output attauiable. 
Hence, any generally applicable and simple expression for the mag- 
netization relations, even if purely empirical, would provide a con- 
venient means for evaluating and describing the performance charac- 
teristics. 

The expressions given below for the magnetization relations are not 
purely empirical, but are those obtained as approximate solutions to the 
magnetic field equations by the magnetic circuit method. When experi- 
mentally evaluated, however, their utility in defining the characteristics 
of the elp(^tromagnet to which they apply is independent of this interpre- 
tation, and depends only on their conformity to the observed magnetiza- 
tion relations. 

It is convenient to formulate the expressions for the magnetization 
curves in terms of the reluctance (R, the ratio JF/^, where fF is the mag- 
netomotive force 4x-V7. The observations of Fig. 1 are plotted in Fig. 3 
in the form of curves giWng (U v.s. .V/ for various values of x. A constant 
value of ip is represented in such a plot by a straight line through the 
origin. The radial lines used as .supplementary co-ordinates are spaced 
to give a convenient scale for (p. 

The reluctance cur\'cs of Fig. 3 are similar in character to those apply- 
ing to most ordinary electromagnets. Each curve has a relatively flat 
characteristic in the ^ncinitv of a minimum located on a common flux 
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line, that marked >p in the figure. At vahies of ^ above this minimum, 
the reluctance increases at an increasing rate, indicating an upper limit 
to the value of ip approached as NI becomes very large. The observed 
value of ip' may be interpreted as that for which the core density cor- 
responds to maximum permeabiUty, while the indicated upper limit may 
be interpreted as the saturation flux ip" . 

The observations plotted in Figs, 1 and 3 were obtained with the 
sample initially demagnetized, a convenient reference condition for 
measurement. In actual use, relays and other electromagnets have been 
previously operated, and the applicable v? versus ff relation is that for 
repeated magnetization. In this case there is little or no increase in re- 
luctance below <p' . For engineering purposes, the reluctance below ip' may 
be taken as constant at the value observed at ip in measurements made 
from an initially demagnetized condition. This assumption is equivalent 
to taking the ^ versus f? relation as linear up to the "knee" of the curve, 
the point of tangency with a line through the origin. 

It follows that expressions for the magnetization relations may take 
(R as a function of % only, independent of ^o, in the low density region, 
^ < tp', while for <p > <p', modified expressions must be employed which 
show (R increasing with ^, and approaching 5/^". 

Magnetic Circuit Schematics 

The reluctance expressions for the magnetization relations can be con- 
veniently defined by the magnetic circuit schematics shown in Figs. 4, 
5, and 6. In these the reluctance (R, or ffV, is represented by a network 
of component reluctances, analogous to a network of electrical resist- 
ances, with flux analogous to current, and magnetic potential analogous 
to voltage. Thus the reluctance corresponding to Fig. 4 is given by: 



(Rl ((Ro + j] 



(R = -^ ^ . (3) 



X 



(Rl + (Ro + J 

The interpretation of the parameters (Ro , (Rz, , and A in terms of the 
magnetic circuit concept is discussed in the companion article cited 
above.^ In the analysis of experimental I'esults they are parameters 
which summarize measurements in which the observed values of (R con- 
form to (3). From this viewpoint Fig. 4 indicates only that the total flux 
<p is the sum of a flux <pi. for which the reluctance is constant, and a flux 
ip(, for which the reluctance varies directly with x. 
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Fig. 4 — Kquivaleut magnetic circuit. 




Fig. 5 — Magnetic circuit for core saturation. 
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Fig. 6 — Miignetic circuit for armature saturation. 
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In Fig. 4, 0lo , (Hi, , and A are constants, so that tfl is a function of x 
only, independent of tp. Equation (3) can then apply only to the low 
density region, tp < tp'. For most ordinary electromagnets, the magnetiza- 
tion relations for <p < ip conform to (3), and hence to the schematic of 

Fig. 4. 

For ip > tp', the expression for (R must provide for the variation with 
<p. It usually suffices to use an expression in which the only term varying 
with <p is that corresponding to the path in which saturation first occurs. 
The most common case is that in which saturation first occurs in the 
core. The reluctance can then be taken as conforming to the schematic 
of Fig. 5, in which mo2 , fllL2 and A2 are constants, and <Rc is a function of 
<p only. The total reluctance is given by: 

(R = (Re + CRb, (4) 

where: 

(JIl2 ( (R02 + -^ 

(H, = — ^- ^ . (5) 

(RlS + {R02 + -J 

If the variation with ^ is taken as conforming to the empirical Frohlich- 
Kenuelly equation/ (ftc is given by: 

(p ~ <p <p — *p 

where tp" is the saturation flux and ip' is the flux for maximum perme- 
ability or minimum reluctance (Re . If these assumptions apply, the 
minimum values of (H for all values of x must lie on ^', as for the results 
of Fig. 3. This common minimum is thus evidence that the core is con- 
trolling with respect to the variation with <p, and that (4) and (0) are 
applicable. .\s (6) applies only for <p > ^', (Re is the vaUie of (Re not 
only at ^ = >p', but throughout the low density region ,p < ip'. In the 
alternative form given by (6), (He is defined by this equation, and repre- 
sents merely the intercept of a plot of (Re extended below the region to 
whi(rh (0) applies. 

In some electromagnets saturation occurs in the armature rather than 
in the core. This is the case, for example, in high speed relays in which 
the armature section is minimized to reduce its mass. In such cases, the 
reluctance conforms to the schematic of Fig. 6(a), in which (R.t represents 
the reluctance of the armature, which varies with (Pa ■ As the ratio ipa/(p 
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then varies witli x, the minimum values of <Si do not lie on a common 
value of <fl in this case. 

In Fig. 6(a), (Rx may be taken as given by an expression of the same 
form as (G), with if' and ip" replaced by the minimum and saturation 
vaKies of ipa , and (Re replaced by the minimum value of (Ra . <Hc may be 
taken either as given by (6), or as a constant if the variation in (R^ is 
dominant. The other parameters of Fig. 0(a) are constants. The expres- 
sion for (R applying to Fig. 0(a) may be obtained in the same maimer as 
those given for the circuits of Figs. 4 and 5. 

The magnetic circuits of Figs. 5 and 6 are called "design" circuits, 
because their con.stants may be estimated from the dimensions and 
material constants of the design, as discussed in the companion article.^ 
Estimates of the corresponding values of the equivalent magnetic circuit 
constants of Fig. 4 are obtained from these design constant estimates by 
the equivalence equations given below. 

Condi/ions of Equivalence 

Whatever magnetic circuit is taken as applying, the component re- 
luctances are independent of tp in the low density region, ip < ip', and 
are therefore constants except for the gap reluctance. If this varies 
directly as x, as assumed in the schematics shown, the reluctance for 
any circuit reduces to an expression of the form of (3), applying to the 
circuit of Fig. 4. As shown in the companion article,^ the reluctance of 
the circuit of Fig. 5 when (Sic is constant is given by (3) when the param- 
eters of this equation arc given by : 

(Rl = (Rt2 + (Re , 

(Ro = P'(R02 + P(R(7 , 

where: 

1 , (Re 
(Rl2 

These relations ai'e the conditions of equivalence, for which the reluc- 
tances of Figs. 4 and 5 are identical for all values of x. With them, the 
reluctance given by (4) can l)e leduced to the simpler form of (3) when 
Sic is a constant, as throughout the low density region. 

Similar relations apply to the magnetic circuit of Fig. 0(a). In particu- 
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lar, when (Ra is a constant, the reluctance (Sip , or Se/<Pa , represented by 
the series parallel path of <pa in Fig- 6(a), may be represented by the 
simple parallel path of ^^ in Fig. 6(b). By analogy with equations (7) 
the constants applying are related by: 

6iLF = (RlA + (Ra , 

A. = ^, (7A) 



where 



flloa = p'(^3 + p(5iA , 



1 _L ^^ 
P = 1 + - 

GIla 



The circuit of Fig. 6(b) is identical mth that of Fig. 5 with 1/(Rl2 equal 
to the sum of I/GIlc and 1/<Rlp ■ Thus the circuit of Fig. 6(a) may be 
reduced to that of Fig. 5 when (Ra is constant. The resulting circuit may 
be reduced in turn to the equivalent circuit of Fig. 4 when (Re is constant. 
Thus the equivalent circuit may be used to represent all cases in the low 
density region, where the component reluctances are independent of ip, 
provided the gap reluctance varies linearly with x. 

3 ANALYSIS OF MAGNETIZATION MEASUREMENTS 

In analysing the observed magnetization relations, it is convenient to 
treat the low density and high density relations in separate and successive 
steps. For the former, the analysis is based on the equivalent magnetic 
circuit of Fig. 4. 

EVALUATION OF EQUIVALENT MAGNETIC CIRCUIT CONSTANTS 

For a given value of x, the total reluctance in the low density region is 
taken as constant at the minimum value 6{' (x) observed in measurements 
made from the demagnetized condition. If the obsei'vations are plotted 
as reluctance curves, as in Fig. 3, these \'alues of dV (x) may be read di- 
rectly. When the recording fluxmeter^ is used, <R' (x) can be obtained 
from the <p versus JF curve as the slope fT,V of the line through the origin 
tangent to the curve. The reciprocals F(x) of the values of (R' (x) thus 
evaluated may be plotted against x. Values of P(x) corresponding to the 
minimum reluctance values of Fig. 3 are plotted In this way in Fig. 7. 
The origin of x is taken as for a gap of 0.025 cm., corresponding to the 
operated position of the actuating card for maximum stop pin height. 
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Values of .t; moasurod from this origin are termed travel, as distinguished 
from gap vahies, which are measured from the position of iron to iron 
contact. 

If the relation between dl' (x) and x conforms to (3), the values of 
P(x) must conform to: 






(8) 



Then if Xc is a central i-efercncc value of .r, and P{xr) the corresponding 
value of P(x), the expression for Pixc) — P(x) given by (8) may be 
written in the form : 



X - Xc 



Thus conformity with (3) requires the ratio given by (9) to vary 
hnearly with x. This ratio is plotted against x — Xc in Fig. 7, referred to 
the upper and right hand scales. The values of this ratio are sensitive to 
deviations in the vahies of x and P(x) near Xc , and minor variations 
from linearity are to be expected here. Allowing for this, the results 
agree with (9), and the relation between (H' (x) and x therefore conforms 
to (3). 

From (9), the slope and intercept of the linear plot may be interpreted 
as indicated in Fig. 7. Then A may be evaluated by dividing the intercept 
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value by the aciuarc of the slope value, autl thi? latter then serves to 
evaluate (ih by .sulitracting Xc/A. On substituting these values of (Ro and 
A in (8) at .t = .r, , oli. may be evaluated. For the case plotted in Fig. 
7, the values thus obtained are: 

(Ro : 0.0300 cm"' 

-4: 1,34 cm' 

iSiz. : 0.0690 cm"' 

These values of the equivalent circuit constants, substituted in (3), 
suffice to accurately characterize the magnetization relations and thus 
the electromagnet's performance through the low density region. 

Evaluation of High Density Relations 

For cases of core saturation, the high density relations may be analysed 
in terms of the magnetic circuit of Fig. 5. The corresponding reluctance 
is given by (4), in which (Re = <Rc for ^ = ^'. If Ql'c can be determined, 
the three parameters dcterminhig (He in (5) may be evaluated by the 
e(iuivalence equations (7). To determine (He through the high density 
region from (0) reeiuires evaluation of (Re, <p', and (p". Thus the analysis 
of the high density relations reciuiros evaluation of these three quantities 
in addition to the three parameters of the e(iuivalent circuit evaluated 
above. 

A method for evaluating (Re , f', and ip" directly from the magnetiza- 
tion curves is described subsequently, following a description of a pre- 
ferred method which requires additional measurements. These measure- 
ments are determinations of the drop ti\c>p in magneti<' potential through 
the core, and are conveniently made with the magnetomotive force 
gauge developed by W. B. EUwood. 

Magnetomotive Force Measurements 

The EUwood mmf gauge measures the difference in magnetic potential 
between the outer ends of the two reeds of a dry reed switch, as deter- 
mined by the coil current required to open the switch. The end of one 
reed is used as a probe, which is brought into contact with the magnetic 
.structure being measured. The end of the other reed, at a distance of 
about 5 cm from the probe, then lies on a relatively distant surface 
surrounding the structure: this surface is substantially at the same po- 
tential for any position of the gauge. Then the difference between two 
such measurements made at different points on the magnetic structure 
measures the difference in magnetic potential between them. 
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Fig. 8 — Relation uf core potciUiiLl drop to applied magnetomotive force. 

To the extent that the physical structure may be identified with the 
corresponding magnetic circuit element, the drop in magnetic potential 
in the core may be identified with Olcv hi Fig. 5. Assuming this identifi- 
cation to apply, magnetic jji'obe measiu-ements may be made at the two 
ends of the core, as at the points marked X in Fig. 2. The potential drop 
observed in these measurements is the magnetic potential ff^ applied to 
the external magnetic circuit. fT/; is equal to the applied magnetomotive 
force 3^(or iirNI) less the potential drop ^c = (Ruip in Ihe core. Thus 
EFc = ff — S?B- Values of tFc/(4T) thus determined for the relay of Fig. 2 
are showai plotted against !T/(47r) for various values of x in Fig. 8. These 
curves are substantially linear in the low density region: their upward 
concavity at higher values of ff is evidence that saturation first occurs in 
the core. 



Evabuidon of Core Relucatance Conslanls 

The magnetization curves of Fig. 1 and the core magnetomotive force 
curves of Fig. 8 were obtained with the same model. For given values of 
J and X there can be read from these two figures corresponding values of 
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if and £f c , giving the forresponding values of core reluctance (Re = ?c/v. 
Such values of J^r have been determined for two values of x (the smallest 
and largest included in these figures), and are shown plotted against 
Jc/(47r) in the dashed curves included in Fig. 3. Values of (Re for inter- 
mediate values of x are intermediate between these two cur^-^es, whose 
similarity supports the assumption of Fig. 5 and equation (4) that (Re 
is substantially independent of x. 

For tp > <p' the plot of (Re versus ffc is substantially linear. In this it 
conforms to (6) which, by substitution of ifc/Slc for tp, may be written 
in the form : 

m, = m^' + ^ = (i-<,V«^ + %- (10) 



tp \ <p / ^ 

Thus ip" may be evaluated from the slope of the plot of (Re versus 
ffc , and (Re from the intercept of the dotted line extension, as indicated 
in Fig. 3. The observed relation deviates from (10) in the vicinity of 
tp' , and the linear relation conforming to (10) does not intersect Olc at 
^', but at some higher value of ^. Thus ((i) more nearly represents the 
observed relation if ^p is taken as this higher value, rather than at the 
flux for minimum (R previously taken as ip' and marked at such in Fig. 
3. The error in using this latter value of ^p' in (G), however, is minor, and 
of little significance in engineering estimates. 

AUemaiive Method of Determining Core Reluctance Constants 

When magnetomotive force measurements are not made, the core 
reluctance constants can be detemiined directly from the magnetization 
curves as follows: 

As can be seen in Fig. 3, the line representing tp" can be directly identi- 
fied as the apparent asymptote of the reluctance curves: specifically by 
determining the flux line parallel to the tangent to the upper portion of 
the lowest reluctance curve. As <p' is the value of ip at which the reluctance 
curves have their common minima, ^' and (p" can be readily evaluated, 
and only (Re remains to be determined. 

On substituting fF/(R for tp in (6), and substituting this expression for 
(Re in (4) the resulting equation can be solved to give: 

(31 = ^(R« + (Re + ^ + /j/^(iU + (Re - P)' + 4(R^ ^) . 

Let (R" be the particular value of (R for which 
g: = (r'^" = («;, + (Rs)/. 
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As both (Rc/(H' and (Rc/Qi are small, substitution of ff 
preceding expression gives as an approximation: 

(ft" = (ft' + Vtftc^t' , 

from which: 

. C(R'' - (ft')' 



(ft ^" in the 



(R' 



(11) 



Thus, as indicated in Fig. 9, the value of (ft' may be read from any 
one of the reluctance curves, and the corresponding vahie of ip"S\.' com- 
puted. The value of <R corresponding to ff = Q{'<p" is read from the curve, 
and taken as (ft". The values of (ft" and (ft' are then substituted in (11) 
to determine (Re- This procedure may be followed for all the reluctance 
curves, and a mean value of (ftc computed, (ftc is then taken as 

/(Rc/(^" - <p'). 

Evaluation of External Reluctance Constants 

The procedures described above pro\nde for the determination (a) of 
the equivalent constants of Fig. 4, (b) of the minimum value (ftc of (R^ 
in Fig. 5. The remaining constants of Fig. 5 can then be computed by 
means of the equivalence equations (7). Table I lists, for the measure- 
ments of Figs. 1 and 8, the values of the equivalent constants deter- 
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Table I — Evaluation of Magnetic Circuit Constants 



Equivalent Values {Fig.4) 


Design Values (Fig. 5) 


From Fig. 7 

(Ho : 0.0300 cin-i 

A: 1.34 cm* 
(Rl : 0.0690 cm-' 


From Fig. 3 

(Re : 0.0070 cm-^ 
From equations (7) 

(Roj : 0.0181 cm-> 
Ai : 1.66 cm^ 

(Rf : 0.0620 cm^i 



mined in Fig. 7, and the value of (Re determined in Fig. 3. These values 
have been substituted in equations (7) to determine the component 
terms of (Rs . 

The constants applying to Fig. 5 are of interest prunarily for compari- 
son \vith the values computed from the design, as discussed in the com- 
panion article.^ While formally required to evaluate the high density 
magnetization relations, they are not explicitly required in estimating 
the high density pull, as is sho^-n in a later section. 

The roluctauce (Re of equations (4) and (5) can be evahiated experi- 
mentaUy through the full range of the magnetization measurements 
when the magnetomotive force measurements are available. Thus Figs, 
i and 3 can be used to deteimine values of ip and ^e (or fF — Jc) for 
corresponding values of fF and x. In this way there have been obtained 
the curves of (Ra plotted against ^e for various values of x in Fig. 10. 
These curves are included here merely to illustrate the approximate 
validity of (4), in which (Re is taken as independent of ip. The departures 
from this condition are minor except for the smaller gaps at high values 
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of If. Here the increase in flux density in the iron parts external to the 
core results in a significant increase in (Sis as tp increases. By comparison 
with the plot of (Sic included in Fig. 1, however, this increase in (Rg is 
minor, and does not affect the fact that the limiting reluctance is 'S/^p" , 
where (p" is the saturation flux of the core. 

It may be noted in passing that the values of (R^ lying on tp must 
conform to (5). As this is of the same form as (3), these values of fllg 
may be used to eA'aluate the component tenns of (o) by the procedure 
applied to the values of (Jl' in Fig. 7. As the same data are emploj'ed, 
the values of fl^02 , ^i.i , and A^ thus determined agree with those com- 
puted by means of equations (7) within the accuracy of the computa- 
tions. 

Case of Armature Saturation 

The analysis of the low density relations is, of course, independent of 
where saturation occurs, and involves merely the determination of the 
equivalent magnetic circuit constants by the procedure described above. 
In the high density region, however, the quantities appearing in Fig. 6(a) 
must be evaluated when armature saturation occurs. This requires 
measurements not only of ^ and <p, but also of 7e and ipA . 

^E is given by the EUwood mmf gauge measurements previously 
described. To measure ^,, re(]uires the use of a search coil wound on the 
armature, located over the region of maximum density. For twin return 
path structures, such as the relay of Fig. 2, twin search coils must be 
used, connected to measure the total armature flux. If the variation of 
(11,1 with <pa is to be determined direcliy, additional mmf gauge measure- 
ments must be made of the potential drop (i{A<PA through the armature. 

With <pa determined as a function of fF^ for various values of x, the 
reluctance (Rp , or ^e '<pa may be analysed by the procedure described 
above for the analysis of the reluctance (R or S/(p. The curves of (Rp versus 
EFb have minima (R'f at tpA . The reciprocals of these values of (Rp are 
plotted and analysed as in Fig. 7 fo evaluate the equivalent constants 
(Sli.p , (Roy , and-l^ of Fig. 0(b). The constants (R1 , ipA , andip" characteriz- 
ing the relation between (ft., and i^.i are determined either from mmf 
gauge measurements or liy the method of Fig. 9. Then the constants 
<5ii.A , <Rti3 , and ,43 of Fig. 6(a) can be evaluated by means of the equiv- 
alence equations (7A). 

To complete the determination of the (]uantities appearing in Fig. 6(a), 
(R/,r is e\'aluated as the ratio ST^ '(<p ~ <pa), while ii\r , which may either 
be constant or conform to (6), is {'\-ahiatcd from the relation between (p 
and EFc . 
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4 PULL RELATIONS 

As shown in Section 2, the pull F can be determined from the field 
energy U by means of eciuation (2). If the gap reluctance varies linearly 
with X, this equation reduces to: 

as shown in the companion article/ Equation (12) is Maxwell's law for 
the pull between parallel plane surfaces of area A. As used here, l/A has 
the more general sense of the coefficient of x in a linear expression for 
the gap reluctance. This equation is therefore applicable to all the cases 
previously discussed, with A, of course, taken as ^2 or A3 for the rela- 
tions applying to Figs. 5 and 6. As with the magnetization relations, it is 
convenient to give separate consideration to the pull at low densities 
and that at high densities. 

Pull for Linear Magnetization 

At low densities, where the magnetization relations are approximately 
Unear, the reluctance confoi-ms to the equivalent magnetic circuit of 
Fig. 4. In this case, ^g is given by ff/(fRfi + x/A), and equation (12) be-, 
comes : 



which may be written in the more familiar form: 

2ir{NIY 

A (m„ + 2) 

For conformity with (13) the pull for a given value of x should vary 
as {Nlf, giving a linear plot with a slope of two when plotted against 
Nl on logarithmic paper. Fig. 11 shows pull measurements of the relay 
of Fig. 2 plotted in this way. The dashed lines tangent to the curves 
have a slope of two, and conform in this respect to (13). It is convenient 
to denote the pull indicated by these dashed lines as F' . Near and below 
the point of tengency the actual pull F differs little from F' . The pull, 
like the magnetization, is measured under the condition of initial de- 
magnetization, which gives lower magnetization, and hence less pull, 
than normally applies in actual use when the magnet has been previously 
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Fig. 11 —Pull curves. 

operated. Hence F' satisfactorily represents the pull applying in practice 
through the low density region. In the high density region, the pull falls 
off as the result of incipient saturation, and the ratio F'/F increases as 
ST (or A-kNI) is increased. 

High Density PuU — Core Saturation 

It is convenient to expres.s the high density pull in terms of the ratio 
F'/F, where /■' is the actual pull, and F' the indicated pull conforming to 
(13) for the same value of ff (or 47rA7). As F'/F is unity at low densi- 
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Fig. 12 — Evaluation of magnetic circuit constante from pull data. ' ' 

Evaluation of Equivalent Magnetic Circuit Constants ' r-y 

The low density pull is given approximately by equation (13), to 
which the dashed F' lines of Fig. 1 1 conform. These lines coincide with 
the actual pull curves at the points of tengency, which must correspond 
to the minimum reluctance values. Taking F as F', (13) may be written 
in the form: 

Values of NI/F' have been detemiined for the dashed lines of Fig. lb 
and are plotted against :t; in Fig. 12 — the plot, marked SF/(47r\/F') ■ T" 
agreement with (18) these points determine a straight line. The slope 
and intercept have Ijecn used to determine the values of (Ro and A listed 
under "Pull Results" in Table II. 

As the reluctance (Rk in the circuit of Fig. 5 is identical in form with 
the reluctance (R of Fig. 4, the pull for Fig. 5 is given by an expression 
similar to (13), with £F, flto , and A replaced by "Se , (R02 , and A^ , respec- 
tively. As the two circuits are equivalent at the points of minimum re- 
luctance and the corresponding points of tangency for the pull curves, 
the pull at these points is the same for the two cases. Thus if the ratio 
of fFs to JF is determined at the points of tangency, a plot of fF.ff/(47r-\/F') 
against x should conform to: 



A^:^/F 



fr = ^^'^ V^'- 



2 V2TA2 



(19) 
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Table II — Evah^ation of Magnetic Circuit Constants 



From Pull Results 


From Magnetization Resujts 


(Ro : 0.0319 cm-' 


0.0300 cm-i 


A: 1 . 40 cm- 


1.34 cm* 


(Ro= : 0.0188 cm-' 


0.0181 cm-' 


Ai : 1.85 cm= 


1.66 rm= 


(Rl : 0.0030 cnr' 


0.0690 cm-' 


(Ri= : 0,0560 cm-i 


0.0620 cm-i 


die : 0.0072 cm-i 


0.0070 cm-' 



The points of taiigency are indicated on the pull curves of Fig. 11. 
These values of 5 have been marked on the corresponding curves of Fig. 
8, from which can he determined the corresponding values of fFc , and 
thus of ^E ■ rt may be noted in passing that these values of Jc are all 
similar, corresponding to a mean level of about 24 ampere turns, or 30 
gilberts. The magnetization results gave values of 4,000 maxwells for 
^' {Fig. 3) and of 0.007 cm"' for (ftc (Table I), giving a value of 28 
gilberts for the product. This agrees with the value of ffc at the tangent 
points of the pull curves, showing that these points coincide with the 
points of minimum reluctance (i^ — ip'). 

Taking the ratio of \9e to fF as read from the curves of Fig. 8 at the 
indicated tangent points, the values of g^/(4ir\/F^) i» Fig. 12 have been 
multiplied by the corresponding values of ^e/^, giving the points indi- 
cated by triangles in this figure. These conform to a straight line relation, 
in agreement with (19). The slope and intercept of this linear plot have 
been used to determine the values of (Ro2 and A2 listed under "Pull 
Results" in Table IT. 

The four quaiititie.s determined from the two plots of Fig. 12 appear 
in the three independent equivalence equations (7), and these three 
equations may therefore be used to determine the other three ([uantities: 
(Rz. , (Ji/.2 , and (iic ■ The resulting values of these latter ([uantities are 
listed under ''Pull Results" in Table II. For comparison, the table in- 
includes the corresponding values of these quantities obtained from the 
magnetization measurements, and gi\-en previously in Table I. The 
agreement between the results derived from these different kinds of 
measurements indicates the validity of the analysis described here. 

High Demn'iij Pvll 

As previously noted, it is convenient to express the pull F in the high 
density region in terms of the ratio F'/F, where F' is the pull computed 
from (13) for the \-aluc of :T applying. F' 'F is readily evaluated from a 
logarithmic plot of the observed pidl versus EF, as illustrated by Fig. 11. 
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Fig. 13 — Pull characteristics for approaching saturation. 



The dashed tangent lines give values of F', so corresponding values of 
F' and F can bp read for various values of 9^. Values of F' /F for the 
results of Fig. 1 1 have been thus determined, and are plotted in Fig. 13 
against (5^/47r) for the values of x applying to the individual curves of 
Fig. 11. 

The vahies of F' /F thus plotted determine a family of straight lines 
having a common intercept at the origin of ff . In this respect these re- 
.sults are representative of the pull characteristics of electromagnets 
when plotted in this way. As an empirical fact, apart from any analysis, 
this linearity provides a useful method of plotting pull observations in 
the high density region, as it facilitates interpolation and the detection 
of inconsistent observations. 

The observed linearity in the relation between F' jF and fF^ conforms 
to the relations postulated in (15) and (17), of which the former should 
apply for core saturation. If (15) applies, the common intercept of 0.65 
at fF^ = in Fig. 13 should eciual 1 — (^'/^") , indicating a value for 
ip' !ip" of 0.59. To determine if the relation fully conforms to (15), the 
relation between the slopes of the lines and the corresponding values of 
X must be compared with that indicated by this equation. 

Let //' bo the ol)served slope of the plot of F' IF versus ff" for a particu- 
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lar value of x. If the relation conforms to (15), ij = 1 '((ilV"). whei-e (ft' 
is given by (3). Writing Cl for the ratio (<Rl + <H(,)/(Ro , and u for the 
ratio x/(A(Ro), 

1 + V 
Cl + U 

and hence the vaUies of y should conform to the equation: 

C, + u 



yd + n) = 



<sW 



(20) 



Using the values of A and (Ro determined from the pull results in Fig. 
12 and tabulated in Table II, values of u have been determined for the 
values of x applying in Fig. 13, and the corresponding values of y have 
been determined from the slopes of the corresponding lines. The values 
of 7/(1 + m) thus detennined are shown plotted against u in Fig. 14. The 
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Fig. 14 — Evaluation of saturation flux v" and leakage reluctance (R^, from pull 
characteristics. 
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Table III — Evaluation of Magnetic Circuit Constants 



From Pull Measurements 


From Magnetiza.tion|^Measurements 


Figs. 13 & 14: 

ip': 4,700 maxwells 
ip": 8,000 maxwells 
(Rl : 0.091 cm-' 
Fig. 12: 

A: 1.46 cm* 
(Ro : 0.0319 cm-' 
CRt : 0.063 cm-' 


Fig. 3: 

tp' : 4,0(H) maxwells 
>p": 8,500 maxwells 

Fig. 7; 

A: 1.34 cm= 
(Ro : 0,0300 cm-' 
(Rl: 0.060 cm^' 



resulting plot is linear, in agreement vnth (20), showing the observed 
relation between F'/F and 5^ to conform fully to (15). 

The expressions fnr the slope and intercept given by (20) are indicated 
in Fig. 14. From the values of the slope and uitercopt given by the plot, 
Cl may be evaluated, giving (Rl , as (Ro is known. With (Rl thus evalu- 
ated, <p" may be determined from the value of the slope. The values of 
(R;, and tp" thus determined from results of Fig. 14 are listed in Table 
III, together with the value of ^' given by the ratio (p'/^p" determined 
from the common intercept in Fig. 13. Included in this table under 
"Pull Results" are the values of A, (Ro , and (R;, determined from the low 
density data and previously tabulated in Table II. Values of the same 
quantities as determined from the magnetization measurements are in- 
cluded for comparison. 

For both pull and magnetization analysis, empirical relations have 
been assumed in (G) and (15) for the variation of SI with tp, agreeing with 
each other and the known conditions only at their upper and lower limits. 
Ill view of this, the agreement in the values of ip' and tp" l)etween the 
pull and magnetization results is good, and the deviation in the Y'alue of 
(Rz, found in Fig. 14 from the other values of this quantity is not surpris- 
ing. 



High Density Pull for Armature Saturation 

The low density pull for armature saturation may be analysed as in 
the case of core saturation. The pull data may be plotted against fF on 
logarithmic paper as in Fig. 11, and the tangents of slope two gi\'ing F' 
dra\vn, as in this figure. The values of ^/-y/Y' may be plotted as in Fig. 
12 to determine .4 and (3lu , and .42 and (R02 also evaluated if "Se has been 
determined from magnetomoti\re force gauge measurements. As before, 
the remaining magnetic circuit constants of Figs. 4 and 5 may be de- 
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termined: the latter in this case are those corresponding to the equivalent 
schematic of Fig. 6(b). 

For the high density pull, values of F'/F may be determined and 
plotted against EF as before. The resulting plot is similar to Fig. 13 in 
showing linear relations for each value of x, radiating from a common 
intercept. For the presentation of pull results, therefore, this method of 
plotting is applicable for both core and armature saturation, and for the 
case where saturation is approached concurrently in both. 

When the slopes if are determined, however, and y(l + u) plotted 
against u as in Fig. 14, a straight line is not usually obtained unless 
saturation is confined to the core. The observed relation for armature 
saturation is usually a curve which is concave upward, and approaches 
the horizontal at small values of u (travel x small). This result can be 
interpreted from the expression (17) for F'/F in the case of saturation 
confined to the armature. The corresponding expression for y is l/ifSiWl)- 
It was shown above that in this case (RJ? is approximately equal to 
mo + x/A or (Ro(l + a). The corresponding expression for y{\ + u) is 
therefore the constant term (l/(mi)^^'). Thus for saturation confined to 
the armature, the plot of y{\ + ") versus u is simply a horizontal line, 
from whose value <pa can be determined. In the more usual case of con- 
current saturation, this relation is approached at small values of u, 
where ipx/^ is relatively large, while at larger values of u, where s^a/v' is 
smaller, the relation approaches that for core saturation. 

6 DISCUSSION 

The procedures for the analysis of pull and magnetization measure- 
-ments described in this article are primarily intended for the evaluation 



of the magnetic circuit constants in development studies. They also have 
some related applications, and the choice of procedure varies with the 
application and with the measurements available. The following discus- 
sion of the apphcations of this analysis indicates the most convenient 
procedures for each case. 

Prosentaiion of Pull Measurements 

Pull measurements are used not only for the guidance of development 
studies, but as engineering data in the application of relays and other 
electromagnets. A convenient form of presentation is that used in Fig. 
11, where the pull is plotted against ampere turns on logarithmic paper. 
In preparing such data, interpolation of the measurements and the 
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recognition of inconsistent observations is facilitated by a supplementary 
plot of F'jF versus £f' (or Q^lf) as in Fig. 13. 

Estimation of Pull 

Estimates of pull are made in preliminary development either from 

estimates of the magnetic circuit constants, computed from the dimen- 
sions as discussed in the companion article, or from magnetization meas- 
urements. The low density pull is given by (13), requiring only the equiv- 
alent circuit constants A and (Ro for its evaluation, and is represented 
by the F' lines of slope two in a logarithmic plot of F versus A^7. Esti- 
mates of the high density pull may be prepared as linear plots of F'/F 
versus (Nlf. As these conform (for core saturation) to (15) they require 
estimates of <p', <p", and of 01', as given by (3). 

Estimates of pull may also be required in preparing puU data for 
engineering use, with respect to the changes in pull associated with di- 
mensional and other variations in a commercial product. Such estimates 
may be obtained by estimating the effect of the variations on the mag- 
netic circuit constants fitting the observed pull, and computing the pull 
for the changed constants in the manner described in the preceding 
paragraph. 

Evaluation of Equivalent Magnetic Circuit Constants 

The equivalent magnetic circuit constants, which define the field 
energy and the pull characteristics through the low density range, 
largely suffice for the prediction of performance. This is illustrated by 
the use of these constants in the analj-'sis of capability relations in one 
companion article^ and in the dynamic relations presented in another. 

These constants are most readily and accurately evaluated from mag- 
netization data by the procedure illustrated in Fig. 7 and the associated 
computations. A check on the values of A and (Ro is provided by analysing 
pull data by the procedure illustrated in Fig. 12. When only pull data 
are available, A and (Ro can be determined in this way, but the deter- 
mination of (Ri requires the measurement of ffg by the Ellwood mmf 
gauge. With CFe known, the pull results suffice for the determination not 
only of (Rl , but of the design circuit constants of Fig. 5 (provided satura- 
tion is confined to the core). 

Evaluation of the Magnetic Design Circuit Constants 

The magnetic design circuit constants are required only for compari- 
son with estimated values, in development studies of the effect on the 
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performance of the configuration, dimensions and materials of the elec- 
tromagnet. They are most readily and accurately determined from coil 
magnetization measurements. It is convenient to use EUwood mmf gauge 
measurements to determine the core reluctance, l)ut these may be 
omitted and the method of Fig. 9 employed, provided saturation is con- 
fined to the core. Subject to this same limitation, the design constants 
may be evaluated from the pull measurements if these are supplemented 
with mmf gauge measurements. 

For armature saturation, the magnetic design circuit constants can 
only be evaluated from, magnetization measurements, which must in- 
clude armature search coil as well as coil flux determinations, and be 
supplemented by mmf gauge measurements. 

Limits of Application 

The validity and usefulness of the procedures described here rest on 
the agreement of measured quantities with the relations used to analyse 
them. As all the relations used lead to linear plots, the extent of agree- 
ment in any specific case is apparent from the plot obtained. So far as 
the presentation of pull results is concerned, this is the only question of 
validity involved. 

The relations found for the magnetization results are used to estimate 
the field energy and the pull. To the extent the magnetization results 
relate to the flux linkages of the coil, the conclusions drawn from rela- 
tions fitting those results are \'alid. This follows from the energy balance 
considerations discussed in Section 2. The supplementary measurements, 
such as those with the mmf gauge or an armature search coil, are in 
principle only convenient means for determining the relations to which 
the coil magnetization conforms. The relations found by these means 
are only valid to the extent of such agreement. 

The conformity of magnetization relations to the expressions used 
here, corresponding to the magnetic circuit schematics, is closest for 
electromagnets in which the reluctance of the iron parts is small com- 
pared with that of the joints and air gaps. This condition is satisfied for 
most ordinary relays and similar electromagnets in the low density 
range, where the expressions given here most closely apply. The expres- 
sions used for the high density range do not give as close agreement, but 
provide a satisfactory basis for engineering analysis, particularly when 
saturation is confined to the core. The treatment is less satisfactory for 
structures that deviate from these conditions, such as those with arma- 
ture saturation, or those with long cores of small cross section, where the 



108 THE BELL SYSTEM TECHNICAL JOXTRNAL, JANUARY 1954 

iron reluctance is a major component through the full range of opera- 
tion. 
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